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Abstract 

* 

Thermal  diffusivity  measurements  have  been  made  using  the  heat  pulse 
method  on  three  3-D  carbon  fibre  carbon  composites , three  1-D  carbon 
fibre  carbon  composites  from  fibres  "F"  and  "A"  and  matrix  graphite. 
Thermal  conductivity  has  been  calculated  from  these  results  and  a 
model  outlined  which  predicts  the  thermal  conductivity  of  a  3-D  composite 
from  the  processed  conductivities  of  the  constituents. 
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1. 


Introduction 


The  thermal  transport  properties  of  three  dimensional  carbon-fibre 
composites  (C.F.C.C.)  have  been  measured  over  the  temperature  range 
300-3000K. 

The  thermal  diffusivity  (a)  has  been  measured  over  this  temperature 

(1  2) 

range  using  the  laser  pulse  technique.  *  Data  thus  obtained  has  been 
converted  to  thermal  conductivity  (X)  using  the  standard  expression 
X  =  a  p  Cp 

where  p  is  density  and  Cp  is  specific  heat. 

The  structure  of  these  materials  is  such  that  it  is  necessary  to 

make  measurements  along  two  axes  in  order  to  characterise  the  properties. 

Although  the  experimental  data  is,  in  its  own  right  extremely 
valuable,  an  additional  aim  of  this  investigation  has  been  an  attempt 
to  model  the  C.F.C.C.  thermal  conductivity  in  terms  of  the  properties  of 
the  constituents.  Three  different  3-D  composites  were  supplied  by  the 
A.F.M.L. ,  Wright-Patterson  A.F.B.  These  were 

a)  Material  A  3-D  orthogonal  C.F.C.C.  using  Fibre 

"F"  material. 

b)  Material  B  3-D  C.F.C.C.  using  Fibre 

"F"  material  using  on  8-harness  satin  weave. 

/ 

c)  Material  C-  3-D  orthogonal  C.F.C.  using  Fibre  "A”  High 

Modulus(H.M. )  fibre  material. 

Samples  of  matrix  material  and  1-dimensional  composites  to  assist  modelling 
were  also  obtained  and  measured  over  the  same  temperature  range. 

Although  the  finished  3-D  C.F.C.C.  consists  only  of  a  porous 
body  containing  carbon  fibre  yarns  and  graphitised  matrix  material  there 
are  many  variables  within  processing  so  an  accurate  description  of  the 
finished  composite  may  be  much  more  complex.  A  fairly  simple  modelling 
approach  has  been  adapted  since  accurately  quantifiable  information  on 
constituent  properties,  from  1-D  composites  and  matrix  material  is  not 


available. 


Data  for  composites,  A  and  B  (fibre  F)  has  been  detailed  in 
the  progress  report  for  the  period  ending  30/9/79.  Consequently  this 
report  will  concern  itself  primarily  with  data  for  material  "C"  and 
comparable  1-D  composites  and  an  evaluation  of  comparisons  with  data 
for  materials  "A"  and  "B". 


2.  Characterisation  of  Materials 

Materials  "A"  and  "C"  are  3-D  orthogonal  C.F.C.C.  in  which  the  carbon 
fibres  are  formed  into  yarns.  In  the  case  of  composite  "A"  this  is 
known  as  a  2-2-3  composite  comprising  two  yarns  of  1440  filament  yarns  in  the 
X  and  Y  direction  and  three  yarns  of  1440  filaments  forming  the  reinforcement 
in  the  Z  direction.  A  unit  cell  of  this  material  is  shown  in  figure  1. 
Composite  "C"  is  similar  but  in  this  case  is  a  2-2-1  composite.  The  X  and 
Y  reinforcement  is  composed  of  two  yarns  of  1000  filaments/site  whereas 
the  Z  axis  reinforcement  is  composed  of  one  3000  filament  yarn.  Material 
B  is  also  a  3-D  composite,  but  in  this  case  the  yarns  are  not  orthogonal. 

The  X  and  Y  yarns  are  woven  into  an  8-harness  satin  weave  to  form  a  rein¬ 
forcement  plane  which  is  pierced  by  the  Z-yarns.  This  type  of 
construction  is  shown  in  figure  2. 

The  manufacture  of  a  typical  composite  begins  with  the  construction 
of  a  fibre  preform.  The  second  stage  is  to  coat  the  preform  with 
pyrolytic  carbon.  This  infiltration  of  the  preform  is  done  using  low 
pressure,  isothermal  deposition  of  vapour  phase  carbon  (C.V.D)  prior  to 
densification  of  the  preform  by  pitch  impregnation,  and  high  pressure 
pyrolysis.  The  C.V.D.  stage  is  instituted  to  stiffen  the  preform  and 
prevent  deformation  during  high  pressure  pyrolysis.  This  impregnation 
process  may  be  typically  repeated  5  times.  The  composite  is  finally 
graphitised  at  3023K  (2750°C).  The  objective  of  this  processing 
is  to  impregnate  the  fibre  bundles  and  fill  crossover  pockets  and  voids 
with  graphitic  matrix  material. 

The  following  additional  specimens  were  provided  by  the  A.F.M.L. 

Wright  Patterson  A.F.B. 

1)  Samples  of  Fibre  "A"  and  Fibre  "F"  material  used 
in  the  composites. 

2)  A  sample  comprised  only  of  graphitised  matrix  material. 

3)  A  1-D  composite  made  using  fibre  "F"  material  that 
had  not  been  subject  to  the  C.V.D.  process  prior  to 
densification  and  graphitisation. 


4) 


A  1-D  composite  made  using  Fibre  "A"  material  that  had  been 


C.V.D.'d  and  graphitised  to  2573K  (2300°C). 

5)  A  1-D  composite  made  using  Fibre  "A"  that  had  been 
graphitised  to  2973K  (2700°C)  but  not  C.V.D.'d. 


Table  1.  Materials  Specification 


3.  Measurement  Technique 

The  U.M.I.S.T.  apparatus  (figure  3)  has  recently  been  described 

( 3 ) 

in  detail  so  will  only  be  briefly  described  here.  The  front  face 
of  a  disc  shaped  sample  is  heated  instantaneously  (<\-  10  3s)  using  a 
100J  ruby  laser.  The  specimen  is  heated  to  the  measurement  temperature 
inside  a  graphite  susceptor  heated  using  an  induction  coil.  The 
subsequent  temperature  transient  on  the  rear  face  of  the  sample  is 
recorded  using  an  infra-red  sensor.  Thermal  diffusivity  is 
calculated  from 


a 


(1) 


where  L  is  the  sample  length. 

t^  is  the  time  taken  to  attain  half  the  maximum  temperature  rise 
w 

—^5*  is  a  constant,  which  in  the  absence  of  heat  losses  =  0.139. 
The  amplified  signal  is  digitally  logged  by  a  microcomputer  which 
samples  data  at  0.001s  intervals. 

Heat  losses  are  inevitable  at  high  temperatures  and  necessitate 

w  (4 ) 

that  — £  decreases.  This  has  been  analysed  by  Cowan  .  To  minimise 

this  correction  and  avoid  overly  large  corrections  for  finite  pulse  time 

affects  it  is  desirable  to  select  a  sample  length  such  that 

0.025s  <  t ,  <  0.2s. Suitable  programming  enables  t  •,  to 'be  computed  to 
5  2 

precision  better  than  10  4s  (including  correction  for  finite  pulse 

time  where  necessary).  The  ratio  of. transient  amplitudes  AT(10t,) 

i 

/AT(ti)  or  AT( 5t  i )  /AT(t-i)  is  then  determined  and  —r  calculated  using 
the  analysis  outlined  by  Cowan.  Finally  the  sample  length  is  corrected 
for  thermal  expansion  and  the  diffusivity  computed  from  these  three 


parameters . 


4. 


Specimens  G.35mm  diameter  and  of  length  chosen  to  fit  the  experimental 
requirements  listed  in  the  preceding  section  were  cored  from  each 
material  billet  provided.  Samples  were  cored  from  the  Z  and  X  directions 
of  the  3-D  C.F.C.C.  materials  since  they  are  symmetrical  in  the  X  and 
Y  directions.  Specimens  perpendicular  and  transverse  were  machined  from 
the  1-D  composites.  For  the  matrix  graphite  the  isotropy  of  the  material 
was  verified  using  samples  cored  from  two  orthogonal  directions.  Sample 
data  is  presented  in  Table  2. 

Density  measurements  were  obtained  using  a  liquid  densitometer 
and  reveal  only  minor  variations  between  samples.  Bulk  densities  are 
in  general  2-3%  lower  which  indicates  a  small  amount  of  open  porosity. 
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Composite 

Type 


3D  CFCC  "A" 


3D  CFCC  "B" 
Fibre  "F" 


i-D  CFCC 
Fibre  "F" 
(2750°C) 


Matrix 

3-D  CFCC  "C" 


1-D  C  FCC 
Fibre  "A" 

CVD  (2300°C) 

1-D  C  FCC 
Fibre  "A  " 

No  CVD 
( 2750°C) 


Table  2.  Specimen  Details 


Orientation  j 

Sample 

No 

i 

.  Length 

1  cm 

Density 
gmcm  3 

r 

X  axis 

1 

0.242 

1.926 

l  2 

0.242 

1.920 

Z  axis 

’  1 

'  0.399 

1.916 

2 

0.398 

1.887 

1 

I 

3 

i  0.298 

1.919 

'  i 

L  .  i 

4 

— 

0.301 

1.925 

1 

X  axis 

0.439 

i  _  •  . 

l  2 

0.439 

9 

~  - - - - 

7  —  — - -  - - 

]  Z  lx  .  s 

1 

0.338 

1.936 

> 

2 

0.336 

1.931 

1 

3 

0.338 

1.928 

L 

r 

Parallel  to 

* 

0.454 

1.910 

fibre  axis 

1  2 

0.458 

1.905 

- - - 

—  -  -  -  -  —  —  -  - 

J  Transverse  to 

1 

0.1995 

fibre  axis 

. 

A 

0.152 

1 

l  ! 

3 

0.1992 

1 

i 

1 

0.297 

1.369 

1 

i 

2 

0.3005 

1.389 

|  X  axis  | 

1 

0.3068 

1.890 

2 

0.257 

1. 876 

— 

.  3  _ 

0.3105 

1.886 

Z  axis 

1 

0.3668 

1.880 

1 

2 

0.3601 

1.888 

Parallel  to  fibre 

X 

0.4039 

1.622 

axis 

2 

0.353 

1.626 

Transverse  to 

1 

0.2065 

1.629 

fibre  axis 

2 

0.2135 

1.624 

Parallel  to 

1 

0.4086 

1.665 

fibre  axis 

2 

0.3561 

1.668 

Transverse  to 

1 

0.2128 

1.701 

fibre  axis 

2 

0.2012 

1.665 

5.  Experimental  Procedure 

Diffusivity  measurements  were  carried  out  on  all  the  materials 
listed  in  section  2  (except  for  the  Fibre  "A"  1-D  composite  graphitised 
to  2573K)  over  the  temperature  range  300-3000K.  It  is  necessary  to  make 
these  measurements  in  four  distinct  stages  as  follows. 

a)  350-750K.  Transients  were  recorded  using  a  Hg^  Cd^_x  Te 
detector  (cut  off  wavelength  5.5  wm)  and  a  synchronous  (detector) 
lock-in  amplifier.  This  was  an  experimental  system  replacing  the 
earlier  InSb  detector  (used  for  fibre  "F"  composites  )  and  was  only 
useable  down  to  approximately  350K  at  acceptable  noise  levels. 

b)  500-1500K.  Transients  were  detected  using  PbS  detector  (cut 
off  wavelength  3.0  pm)  and  the  measurement  temperature  recorded  using 
a  Type  K  thermocouple. 

c)  1200-2300K.  As  b)  above  but  temperatures  recorded  using  an 
optical  pyrometer. 

d)  2000-3000K.  A  vacuum  (<  10  4  Torr)  was  sufficient  to  protect  the 

samples  over  the  preceding  temperature  ranges.  Above  2300K  the  vapour 
pressure  of  graphite  increases  so  it  is  necessary  to  suppress  evaporation 
This  was  done  by  completing  measurements  in  an  atmosphere  of  helium 

at  20  p.s.i.  pressure. 

The  above  temperature  ranges  all  overlap.  This  is  to  confirm  that 
these  changes  in  measurement  technique  did  not  yield  different 
diffusivity  values. 

For  modelling  purposes  it  is  necessary  to  convert  the  measured 
thermal  diffusivity  to  thermal  conductivity.  This  requires  accurate 
data  for  specific  heat  and  thermal  expansion.  The  specific  heat  data 

t  C  \ 

used  was  supplied  by  A.F.M.L.  synthesised  from  measurements  from 

( 8 ) 

...  (7)  and  a  C.F.C.C.  composite  "A"  from 

350-1000K  on  P0C0  graphite  r 

1500-3000K.  This  data  is  fitted  to  a  5th  order  polynomial  using  a 
least  squares  function  (figure  3). 


Cp  =  -5.444xlO_1  +  5.5076xl0~3T  -  4.9454xl0_6  T2 

+2.3389xlO_9T3  +  5 . 5749xlO~ 1 3T4  -  5. 3241xlO~ 1 7T5 

(  0  ) 

Thermal  expansion  data  supplied  by  A.F.M.L.  has  been  used  to 
correct  for  density  changes  with  temperature. 

Specimens  for  microstructural  investigation  were  prepared  by 
grinding  and  polishing  down  to  1/4  pm.  These  were  then  etched  for 
1-5  minutes  in  a  hot  solution  of  potassium  dichromate  in  phosphoric 
acid.  In  order  to  improve  resolution  in  the  electron  microscope, 
the  etched  specimens  were  then  lightly  sputtered  with  gold . 
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Results 


At  least  two  specimens  were  measured  for  each  type  of  specimen. 
Conductivity  values  derived  from  the  diffusivity  data  were  fitted 
with  least  squares  functions.  A  number  of  different  functions 
were  tried  and  the  best  fit  over  the  whole  temperature  range  was 
usually  a  fourth  order  polynomial.  Although  data  for  each  specimen  type 
is  presented  briefly  in  section  6(A)  -  6(G)  the  coeficients  and  root  mean 
square  errors  for  the  best  fit  functions  are  listed  in  Table  3. 

(A)  3-D  Composite  'A* 

(i)  X  Axis 

The  diffusivity  of  two  X  axis  samples  was  measured  over  the  complete 

temperature  range.  The  results  are  shown  in  Figure  4.  They  provide  an 

excellent  agreement  with  AFML  flash  diffusivity  data,  also  shown,  which 

was  obtained  from  a  CFCC  similar  to  ’A'. 

Figure  5  shows  the  X  axis  conductivity  together  with  directly 

measured  data  obtained  from  different  billets  of  'A'  type  material  by 

(  6  ) 

the  comparative  rod  method'  .  The  direct  method  conductivity  is 
higher  at  low  temperatures,  and  this  may  simply  be  a  product  of  inter¬ 
billet  variations.  However  Lee  and  Taylor^^  and  Minges^10^  have 
observed  similar  discrepancies  between  direct  and  diffusivity-derived 
conductivity. 

There  is  no  systematic  difference  in  the  properties  of  the  two 
samples  and  this  is  shown  more  clearly  in  Figure  6  which  is  a  plot  of 
%  deviation  from  the  least  squares  function. 

(ii)  Z  Axis 

Figure  7  shows  the  diffusivity  of  four  samples,  again  compared 
with  diffusivity  data  from  a  material  similar  to  composite  ’A'^&^. 

As  in  the  previous  case,  the  agreement  of  the  diffusivity  data  is  good. 

Figure  8  shows  the  resultant  conductivity  values  with  comparative 
rod  data  for  comparison^.  As  before,  the  latter  data  is  higher.  Billet 

(9) 

variations  or,  as  Lee  and  Taylor  imply,  differences  in  measurement 


method  may  be  responsible.  In  this  case,  however,  it  seems  appropriate 
to  question  the  validity  of  the  low  temperature  specific  heat  data.  Over 
a  temperature  range  of  350-1000K  this  relies  entirely  on  the  specific 
heat  of  POCO  graphite,  whereas  there  is  evidence  that  the  specific  heat 
of  fibre  'F'  alone  is  30%  higher  at  300K  .  The  law  of  mixtures  would 

suggest  that  the  composite  value  should  lie  between  these  two  limits. 

Figure  9  shows  the  data  deviation  plot  for  the  least  squares  function 
given  in  Table  3.  A  systematic  relationship  between  deviation  and  sample 
can  be  seen  which  may  be  partly,  but  not  completely,  explained  by 
differences  in  sample  density  (see  Table  2). 

(B)  3-D  Composite  'B' 

(i)  X  Axis 

Diffusivity  data  and  derived  conductivity  values  of  two  X  axis  samples 
are  shown  in  Figures  10  and  11  respectively.  No  comparative  diffusivity 
data  is  available.  However,  comparative  rod  conductivity  data^  , 
from  a  different  billet,  shows  excellent  agreement. 

Figure  12  shows  the  deviation  from  the  least  squares  fit. 

(ii)  Z  Axis 

Diffusivity  and  conductivity  results  from  three  samples  are  shown 
in  Figures  13  and  14  respectively.  Comparative  rod  data^&\  again 
from  a  different  billet,  is  higher  by  a  factor  of  20-25%  over  most  of  the 
temperature  range. 

Figure  15  shows  that  deviation  from  the  least  squares  fit  is 
broadly  sample  dependent.  Differences  in  sample  densities  are  insignificant. 

(C)  3-D  Composite  1 C 1 
(i)  X-Axis 

Diffusivity  data  and  derived  conductivity  values  for  three  X-axis 
samples  are  shown  in  figures  16  and  17.  No  comparative  data  is  available. 
There  was  some  indication  of  difference,  between  sample  1  and  samples 
2  and  3.  This  data  is  very  similar  to  that  obtained  for  composite  ’A' . 


Figure  18  shows  the  deviation  from  the  least  squares  fit  given 
in  Table  3. 

(ii)  Z  Axis 

Diffusivity  data  and  derived  conductivity  values  for  two 
Z-axis  samples  are  shown  in  figures  19  and  20  respectively.  There  is 
no  detectable  sample  difference  and  again  similarity  to  the  data  obtained 
for  composite  'A'. 

Figure  21  shows  the  deviation  from  the  least  squares  fit  given 
in  Table  3. 

(D)  Matrix  Material 

Diffusivity  and  conductivity  results  from  two  bulk  matrix  samples 

are  shown  in  Figures  22  and  23  respectively. 

The  samples  were  cut  orthogonally  to  determine  the  degree  of 

anostropy,  if  any.  The  data  indicate  an  essentially  isotropic  material. 

The  high  porosity  of  the  bulk  matrix  (^40%)  contrasts  with  the  high 

density  of  the  composites  and,  by  inference,  with  the  low  porosity  of  the 

matrix  material  within  the  composite.  Evidence  suggests  in  general 

that  the  density  of  carbon  composite  phases  in  situ  is  higher  than  when 

processed  individually^'*''*'^.  However,  the  matrix  ’crossover’  pockets  (see 

Figure  1)  are  known  to  have  porosity  values  of  between  10-50%^6\ 

The  matrix  conductivity  is  very  similar  to  that  of  P0C0  graphite 

type  AXM-5Q1,  as  well  documented  reference  material^  .  With  an  average 

density  of  1.75  grms/cm3,  the  P0C0  poro’sity  is  considerably  lower, 

however.  Whilst  radiative  and  gaseous  conduction  modes  may  be.come  important 

(12 ) 

in  high  porosity,  low  conductivity  materials  no  enhancement  of  the  matrix 

thermal  properties  was  observed  in  a  helium  atmosphere  both  at  low  and  high 
temperatures.  Similarly,  such  expressions  as  are  available  to  quantify 
raidative  conduction  in  porous  materials  suggest  that  this  too  would  be 


insignificant  in  the  present  case. 


(E)  Fibre  "F”  1-D  Composite.  No  CVD  -2700  C  Graphitisation 


(i) 


Axis  Parallel  with  Fibres 

Diffusivity  and  conductivity  results  from  two  parallel  (//)  axis 


samples  are  shown  in  Figures  24  and  25.  No  comparative  data  is  available. 

The  1-D  composite  was  unique  in  that  its  properties  change 
significantly  after  heating  to  temperatures  >  2300K.  Spot  measurements, 
at  temperatures  of  600-700K,  were  made  after  the  principal  measurements 
had  been  completed.  These  revealed  that  the  diffusivity  of  both  samples 
was  reduced  to  50%  or  less  of  its  original  value.  In  addition  there 
was  a  permanent  change  in  simple  length  and  diameter  of  +  1-1.5%  respectively. 

Since  there  are  no  orthogonal  yarn  influences,  it  is  expected 
that  the  1-D  composite  thermal  expansion  be  higher  than  that  of  the 
3-D  materials'  .  However,  no  permanent  offsets  have  been  reported 
before.  Significantly,  it  is  now  known  that  the  processing  of  the 

1- D  material  differs  from  that  of  the  other  composites''  reported  here. 

Typically,  the  high  temperature  (2000-3000K)  measurement  runs  lasted 

2- 3  hours.  Oxygen  contamination  has  been  ruled  out  since  other  sample  types 
were  not  effected. 

Subsequent  liquid  densitometer  measurements  showed  a  density  increase 
of  about  1%.  Compared  with  the  nominal  decrease  in  bulk  density,  this 
indicates  a  decrease  in  closed  porosity  and  micrstructural  changes. 

The  high  conductivity  of  the  1-D  parallel  axis  is  another  surprising 
feature.  The  room  temperature  value  of  approximately  3.0  W/cmK  must  be 
compared  with  available  data  on  Fibre  ’F'  conductivity  which  indicates  a 
value  of  =  0.6  W/cmK  along  the  fibre  axis^*^’^  .  This  large  difference 
in  the  properties  of  the  1-D  composite  .and  fibres  has  important 
implications  and  will  be  discussed  more  fully  when  composite  modelling 
techniques  are  examined. 

Doubts  have  been  raised  concerning  the  efficacy  of  diffusivity 
measurements  on  highly  orientated  composites  where  the  integrity  of 
conduction  paths  is  preserved.  A  good  example  of  this  is  the  parallel 
axis  of  the  1-D  composite.  Arguably,  an  effective  composite  diffusivity 
does  not  exist  for  such  materials  in  which  event  the  rear  face  temperature 
analysis  (see  Equation  2)  will  give  meaningless  results.,  except  in 


two  extreme  cases. 
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These  are: 

u2 
a2 

where  and  are  the  diffusivities  of  the  phases. 

(9) 

Case  2)  has  been  used  by  Lee  and  Taylor  to  measure  the  diffusivity 
of  carbon  fibres. 

The  rear  face  temperature  analysis  enables  the  calculation  of 
diffusivity  with  any  value  of  x  (0<x<l)  in  t^provided  the  appropriate  value 
of  w/tt2  is  used  in  Equation  2.  Calculation  of  diffusivity  for  different 
values  of  x  gives  a  check  on  how  closely  the  observed  transient  conforms 
to  its  theoretical  form  and  thus  may  provide  information  as  to  the  nature 
of  the  sample. 

The  rear  face  temperature  transients  of  several  1-D  composite  samples 
were  recorded  over  a  temperature  range  where  heat  losses  would  not  distort 
the  transient.  The  difference  in  diffusivity  values  calculated  from  t^ 
values  in  the  range  0.2  <  t^  <  0.8  was  less  than  ±4%.  It  may  be  concluded 
therefore  that  the  parallel  axis  does  have  a  meaningful  effective 
diffusivity. 

(ii)  Axis  Transverse  to  Fibres 

Diffusivity  and  conductivity  results  for  the  transverse  fibre  (TF) 
axis  from  three  samples  are  shown  in  Figures  26  and  27. 

The  samples  transverse  to  the  fibre  axis  also  showed  irreversible 
dimensional  changes  after  heating  to  2400K.  Increases  in  sample  length 
were  higher  being  in  the  range  2. 6-4. 9%.  Fibre  thermal  expansion  is 
known  to  be  higher  in  the  transverse  direction^’ ^ . 

(F)  1-D  Composite  Fibre  "A".  No  CVD  -  Graphitised  at  2750°C 

(i)  Axis  parallel  to  fibres 

Measured  thermal  diffusivity  and  derived  thermal  conductivity 
data  for  two  samples  are  shown  in  figures  28  and  29.  The  diffusivity 
is  some  30%  lower  than  the  data  obtained  for  the  fibre  "F"  -  1-D  composite. 

The  density  is  also  significantly  lower  (1.7  gm  cm  3  compared  to  1.89  gm  cm  3  ) 


1)  a. 


>> 
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and  the  derived  conductivity  values  lie  some  37%  lower. 

(ii)  Axis  transverse  to  fibres 

In  similar  fashion  the  thermal  diffusivity  data  obtained  for  this 
direction  (figure  30)  are  also  some  30%  lower  and  the  thermal  conductivity 
values  (figure  31)  are  likewise  37%  lower. 

(G)  1-D  Composite  Fibre  "A11  +  CVD  Graphitised  at  2300°C 

(i)  Axis  parallel  to  fibres 

Thermal  diffusivity  data  for  2  samples  are  shown  in  figure  32. 

These  data  are  much  lower  than  those  for  the  1-D  composite  of  Fibre  "A" 
without  CVD  by  a  scaling  factor  of  0.3-0. 5.  The  influence  of  CVD  is  almost 
impossible  to  ascertain  since  the  difference  is  almost  certainly  due 
primarily  to  the  difference  in  graphitisation  temperature.  An  attempt  was 
made  to  graphitic  one  sample  by  heating  to  3000K.  However  the  diffusivity 
decreased  and  subsequent  examination  revealed  cracks  between  fibres 
and  between  fibres  and  matrix. 

(ii)  Axis  transverse  to  fibres 

Thermal  diffusivity  data  are  presented  in  figure  33  obtained  for 
this  direction  showed  values  exhibiting  a  different  temperature  dependence 
to  the  values  obtained  for  fibre  "A"  with  no  CVD.  At  T  >  1000K  values 
are  similar  but  at  300K  values  are  only  70%  of  the  non  CVD'd  samples. 


lft. 


7 .  Micron  lructur.il  Invest  if, at  ion. 

Figures  34  and  35  show  scanning  electron  microscope  (C.E.M.)  views 
across  the  fibre  axis  of  two  1-D  composite  samples,  made  from  fibre  "F" 
material.  The  latter  is  sample  Z  after  measurements  up  to  2900K  whereas 
the  former  lias  undergone  no  heating.  There  were  no  obvious  microstructural 
changes  which  could  explain  the  large  transverse  off sett.  Fibre/matrix 
interfaces  appeared  very  similar  with  no  evidence  of  increased  separation, 
as  higher  magnification  shows  (figures  36  and  37). 

Figure  38  shows  a  general  view  across  the  X  axis  of  an  unheated 
specimen  of  C.F.C.C.  "A".  The  fibre  yarns  running  from  bottom  left  to 
top  right  are  Y  axis  reinforcement,  the  other  fibre  yarns  are  Z  axis.  Pores 
of  up  to  lOOwm  were  visible  in  the  matrix  "cross  over"  pockets. 

Broken  yarn/yarn  interfaces  were  a  common  feature  of  the  3-D  composites 
as  illustrated  in  figure  39  whereas  yarn/matrix  interfaces  were  generally 
intact .  The  X  axis  reinforcement  is  composed  of  two  yarns  but  there  was 
little  evidence  of  yarn  splitting  in  the  unheated  composite. 

Figure  40  shows  greater  separation  between  fibres  and  matrix  than 
was  visible  in  the  1-D  composite.  In  some  regions  the  volume  between 
adjacent  fibres  was  completely  filled  with  material  of  an  apparently  different 

form  to  the  rest  of  the  matrix  (see  figure  41).  This  quite  possibly  is  CVD 

+  .  .(6) 

material 

A  general  view  across  the  Z  axis  of  a  C.F.C.C.  "A"  specimen  that 
had  been  measured  up  to  2900K  is  shown  in  figure  42.  Splitting  within  the 
Z  axis  yarns  could  be  seen  and  the  pattern  of  broken  yarn/yarn  interfaces 
repeated  but  more  clearly  marked  than  in  the  unheated  specimen  (figure  43). 
However  spot  diffusivity  measurements  taken  during  cooling  had  revealed  no 
measurable  change  in  thermal  properties.  No  evidence  could  be  seen  of 
changes  in  the  fibre /matrix  interface  (figure  44). 

A  general  view  across  the  X  axis  of  C.F.C.C  .  "B"  is  shown  in  figure 


45,  illustrating  the  different  method  of  fabricating  this  composite.  The 
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X-reinforcements  were  not  easy  to  distinguish  individually  because  of 
the  apparent  absence  of  matrix  "cross  over"  pockets.  The  area  between 
the  Z  yarns  appeared  very  homogeneous  a  feature  further  illustrate  d 
in  figure  46.  Presumably  this  is  a  result  of  using  a  weave  instead  of 
orthogonal  fibres  held  in  place  before  CVD  by  a  jig^\  The  X  and  Y 
reinforcements  are  therefore  able  to  expand  into  the  cross  over  pockets, 
otherwise  filled  by  the  matrix.  The  vertical  cracks,  both  major  and  minor 
were  regularly  spaced  and  probably  mark  the  boundary  between  adjacent 
X-Y  reinforcement  planes  which  are  nominally  0.254mm  (254pm  thick). 

Broken  Z  yarn/X-Y  interfaces  were  quite  prominent  and  the  fibre/matrix 
microstructure  was  very  similar  to  that  of  C.F.C.C.  "A".  Figure  47  shows 
a  general  view  across  the  Z  axis  illustrating  the  fine  weave  pierced  fabric 
construction.  Broken  yarn/yarn  interfaces  are  clearly  visible. 

However  when  we  examine  the  composites  nominally  made  from  fibre  "A" 
then  discrepancies  become  evident.  Sane  very  good  micrographs  were  obtained 
of  the  two  1-D  fibre  "A"  composites  which  clearly  show  a  circular  fibre 
shape.  In  figure  48  for  the  1-D  composite  that  has  been  CVD’d  can  be  distinctly 
seen  a  number  of  spherical  fibres  of  about  6pm  diameter.  At  higher 
magnification  in  a  more  heavily  etched  region  of  the  specimen  (figure  49)  can 
be  seen  two  fibres  which  have  clearly  visible  around  them  a  sheath, 
presumably  of  CVD  graphite.  In  figures  50,  51  and  52  at  magnifications  of 
x  2,000,  x  5,000  and  10,000  are  three  more  regions  of  the  CVD'd  1-D 
composite.  Again  the  graphite  sheath  around  each  individual  fibre  is 
evident  but  more  significantly  the  structure  of  the  in-fill  graphite,  presumably 
matrix  appears  to  orient  itself  transverse  (TOG)  with  respect  to  the  sheath 
of  CVD  graphite  around  each  fibre  which  appears  to  by  parallel  oriented  graphite 
(P.O.G. ) 

By  way  of  contrast  the  non  CVD  unidirectional  composite  has  a  markedly 
different  microstructure.  In  figure  53  is  shown  a  micrograph  at  x  2,000 
which  may  be  compared  with  figure  50.  Again  whilst  the  cylindrical  nature  of 
the  fibres  is  clearly  evident  there  is  a  relative  lack  of  detail  regarding 


the  microstructure  of  the  matrix.  Progression  to  higher  magnification 


figure  54  (x  5,000)  and  figure  55  (x  10,000)  which  are  to  be  compared  with  ' 

i 

figures  51  and  52,  reveal  little  more  detail.  At  x  20,000  magnification 
(figure  5G)  some  detail  is  revealed  with  slight  evidence  that  matrix  j 

graphite  does  form  a  sheath  around  the  fibres.  However  no  further  detail  ' 

can  be  drawn  regarding  the  nature  of  the  microstructure.  When  a 
longitudinal  section  of  the  non  CVD'd  1-D  fibre  "A"  composite  is  examined 
a  little  more  detail  is  apparent.  In  figure  57  there  is  clear  evidence  that 
the  matrix  graphite  surrounding  the  fibre  is  transversely  oriented. 

If  we  now  consider  the  3-D  composite  supplied  as  being  made  of 
fibre  "A"  material  then  a  disconcerting  difference  becomes  apparent.  In 
figures  58  and  59  are  electron  micrographs  of  the  Z  axis  direction  and  X 
axis  direction  respectively.  Higher  magnifications  however  clearly  show 
that  the  fibres  are  crenulated  (figure  60).  This  specimen  has  been  heavily 
etched  to  show  the  relief.  Even  on  more  lightly  etched  samples  figures 
61  (x  5,000)  and  62  (  x  10,000)  the  fibre  shape  can  clearly  be  seen.  In 

these  two  specimens  can  also  be  seen  some  outcrops  of  matrix  graphite 
which  has  obvious  structural  features  of  similarity  with  the  CVD'd  1-D 
composite  (figure  50). 

Figures  63  (  x  200)  and  64  (  x  1,000)  exhibit  evidence  of  cracks 
at  yarn/yarn  interfaces  whereas  there  is  a  better  continuity  between  the 
interfaces  between  matrix  and  yarns.  There  is  evidence  of  substantial 
porosity  in  the  matrix  graphite  in  figure  64. 

It  is  clear  that  the  3-D  composite  is  not  made  from  fibre  "A" 
unless  some  extremely  unusual  etching  effect  is  evident  in  the  3-D  composite 
that  is  not  apparent  in  the  1-D  composite.  We  must  therefore  conclude 
that  the  wrong  composite  has  been  supplied  and  moreover  is  most  probably  made 
from  fibre  "F"  material.  This  does  not  appear  to  be  composite  "A",  although 
the  diffusivity  results  are  very  similar.  Analysis  of  figures  58  and  59  suggests 
that  the  unit  cell  dimensions  in  the  Z  and  X  directions  are  0.685mm  and 

590mm  respectively  which  is  substantially  less  than  the  unit  cell  dimensions 

A 


of  composite  "A."  (0.8't  x  0.76  mm).  This  is  a  scaling  factor  of  0.8 
and  it  is  extremely  unlikely  that  the  electron  microscope  calibration 
could  be  in  error  by  20%. 


8 .  Modelling  of  C.F.C.C.  Thermal  Conducitivty 

The  intended  approach  of  deterrning  the  thermal  transport 

properties  of  the  individual  fibre  and  matrix  phases  and  combining  this 

data  in  the  appropriate  3-D  geometry  is  clearly  invalidated  by  the  lack 

of  suitable  specimens  from  which  1-D  composites  and  matrix  graphite 

experimental  data  may  be  analysed,  and  by  the  fact  that  composite  "C"  does 
not  appear  to  be  made  from  fibre  "A"  material. 

The  thermal  conductivity  of  a  composite  will  be  a  function  of 

the  thermal  conductivities  of  the  matrix  material  and  of  the  fibres. 

These  will  be  affected  by  processing  variables;  matrix  graphite  by 
the  presence  of  porosity  and  the  fibres  by  preform  stiffening  by  carbon 
vapour  deposition.  To  that  extent  none  of  the  back-up  materials;  the 
three  1-D  composites  and  the  matrix  graphite,  can  be  said  to  be  represent¬ 
ative  of  the  3-D  composites. 

The  bulk  density  of  the  matrix  graphite  is  only  1.36  gmem  3 

whereas  the  bulk  densities  of  the  composities  are  much  higher.  This 

inplies  that  some  scaling  factor  is  required  to  predict  the  conductivity  of 

the  matrix  material  of  the  composites.  Again  whilst  composites  A  and  B  and 

the  1-D  fibre  "F"  composite  have  comparable  densities,  that  for  composite 

C  and  the  two  fibre  "A"  1-D  composites  differ  significantly.  Additionally 
the  fibre  F  and  fibre  A  composites  graphitised  to  2750°C  did  not  have 

CVD  treatments  and  cannot  be  consider  typical  of  the  3-D  composites 

whereas  the  only  1-D  composite  that  did  have  CVD  treatment  was  only 

graphitised  to  2300°C.  Hence  it  is  not  possible  to  take  into  account 
the  influence  of  processing  upon  the  properties  of  the  constituent  phases. 

This  may  be  qualitatively  summarised  as  follows: 

1)  Constituent  densities  tend  to  be  significantly  higher  within  the 

CFCC  Increases  of  20%  and  more  have  been  reported  in  both  matrix 

and  fibres^^.  In  the  case  of  the  fibres,  some  of  this  density  increase 

(14) 

is  due  to  the  filling  of  pores  in  the  fibre  bundles 

2)  If  used,  the  CVD  process  produces  a  sheath-like  coating  around  the 
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individual  fibres  (see  Figure  31).  The  form  of  the  coating  varies  with 


fibre  type  and  in  the  case  of  fibre  ' F'  is  isotropic  or 


transverse  .  Where  fibres  are  closely  packed  (separation  less  than  l\ 


fibre  diameters)  the  inter-fibre  space  may  be  entirely  filled  with  CVD 

.  ,  (17) 
material 


3)  The  elastic  modulus,  thermal  conductivity  and  bulk  density  of  CVD 
CFCC  are  lower  than  when  CVD  is  absent.  The  CVD  process  tends  to  seal 


off  the  fine  porosity  of  the  fibres,  preventing  infiltration  by  the  matrix 
(14,15,16) 


4)  The  matrix  crystallites  in  the  vicinity  of  fibres  with  isotropic  CVD 

(14) 

coatings  are  transversely  oriented  (TOG)  w.r.t.  the  fibre  axes 

In  other  regions  e.g.  cross-over  pockets  and  where  spacing  between  fibres 

is  large,  the  matrix  is  isotropic. 

5)  Matrix  crystallites  in  the  vicinity  of  non  CVD'd  ' F*  fibres  are 

(’7) 

parallel  oriented  (POG)  with  the  fibre  axes 

It  is  clear  that  the  properties  of  CFCC  materials  result  from 
a  complex  interaction  of  the  individual  phases.  In  the  case  of  non 
CVD'd  yarns,  the  matrix  must  not  only  increase  fibre  conductivity  by 
infill  of  fibre  defects,  but  also  with  its  POG  structure  contribute 
very  significantly  to  axial  conductivity. 

In  the  case  of  CVD'd  material,  it  is  more  likely  that  CVD  infill 
is  the  major  influence  on  axial  conductivity  with  that  of  the  matrix  less 
prominent.  In  all  examples  however  the  matrix  properties  themselves  will 
be  influenced  by  FVF  since  this  is  a  factor  in  determining  the 
proportions  of  POG/TOG  and  isotropic  matrix.  A  schematic  model  of  composite 
A  is  reproduced  in  figure  65. 
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8.1  Matrix  Graphite 

The  large  difference  in  density  between  the  as  supplied  matrix 
graphite  and  the  fibre  composites  raises  doubts  as  to  whether  the 
measured  values  are  appropriate  for  use  in  composite  modelling. 

It  has  been  conclusively  shown  that  the  conductivity  of  graphite  is 
dominated  by  conduction  along  layer  planes  and  that  the  thermal 
conductivity  of  polycrystalline  graphite  may  be  represented  by  '  ’ 


i  =  £S  r  =  »  r  <2> 

a  a 

where  X^  is  the  thermal  conductivity  along  layer  planes  and  the 

parameters  e  and  8  represent  a  tortuosity  factor  and  a  porosity  factor 

(18) 

respectively.  Taylor  et  al  have  analysed  a  number  of  graphites  using 
this  form  of  analysis  and  consider  that  three  scattering  processes 
contribute  to  the  thermal  resistivity. 


where  X^ ,  X^,  and  X^  represent  the  contribution  to  thermal  conductivity 

due  to  Umklepp  scattering,  grain  boundary  and  isotope  scattering. 

Taylor  et  al  have  shown  that  that  isotope  scattering  will  contribute 

some  2-3%  to  the  total  thermal  resistance  so,  in  view  of  the  assumptions 

to  be  made,  this  will  be  neglected.  An  empirical  relationship  for  the  mean 

(19) 

free  path  has  been  derived  by  Taylor 

Xu  =  8.75xl0"7  exp  1130/T  *  (4) 


from  which  the  thermal  conductivity  Xu  =  1/3  CVL>u  may  be  calculated . 
(Values  are  listed  up  to  1000K  by  Taylor).  Explicit  equations  for 
boundary  scattering  have  been  derived  by  Kelly 


Kl  + 


K2  + 


(5) 


T 


+ 


T 
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whex-e  the  symbols  No,V  ,V  ,h,6,o  ,0. ,0  and  W3  are  as  used  previously 

X  cl  la  T 

by  Kelly. 

At  high  temperatures  it  is  reasonable  to  assume  that  UmkJapp 
scattering  will  dominate  and  that  boundary  scattering  will  contribute 
<  10%  to  the  total  thermal  conductivity.  Data  for  matrix  graphite 
can  be  fitted  by  a  curve  for  which  a  =  9.0  and  =  2,750  °A. 
to  an  accuracy  of  ±  8%  over  all  the  curve.  This  value  of  crystallite 
size  is  slightly  higher  than  those  derived  by  Taylor  for  a  range  of 
polycrstalline  graphites  but  similar  to  a  value  determined  by  Kelly 
and  Gilchrist^ 21^  for  pyrolytic  graphite  graphitised  at  2750°C.  The 
value  of  9.0  for  the  porosity /tortuosity  faccor  is  higher  than  those 
•noted  by  Taylor;  however  the  low  density  (40%  porosity)  of  the  bulk 
matrix  graphite  is  probably  a  prime  contributory  factor. 

This  matrix  graphite  with  a  bulk  density  of  1.36  gm  cm  3  has  to  be 
compared  with  bulk  densities  in  the  range  1.88-1.92  gm  cm-3  for 
the  two  1-dimensional  fibre  "A"  composites.  Hence  to  use  the  raw 

derived  thermal  conductivity  data  is  probably  unrealistic.  A  fully 

dense  isotropic  graphite  would  have  zero  porosity  factor  and  a 

tortuosity  factor  of  /2~  .  A  line  of  interpolation  between  a  =  1.41 

at  zero  porosity  and  a  =  9.0  at  40%  porosity  should  however  permit  an 

evaluation  of  a  for  any  value  of  density.  Although  there  is  'no 

justification  for  this, it  is  perhaps  a  reasonable  assumption  in  view  of  the 

complex  interaction  of  variables  and  lack  of  any  available  data. 
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8.2.  1-D  Composites 

The  thermal  conductivity  of  a  unidiroctionally  rcinfored  composite 
in  a  direction  parallel  to  the  fibres  X^  1  is  usually  expressed  in  terms 
of  the  volume  weighted  conductivities  of  the  constituents  using  an  Ohm's 
law  approach 

X  11  =  V,  X  11  +  (1-V_)  X  (8) 

c  r  f  f  m 

where  X^-11  is  the  conductivity  of  the  fibres 

X  is  the  matrix  conductivity 
m 

is  the  fibre  volume  fraction 

Applying  the  results  for  the  three  unidirectional  composites  should  in 

principle  permit  the  calculation  of  the  thermal  conductivity  X-11  of  the 

fibres.  However  the  significantly  higher  densities  of  the  1-D  composites 

suggests  that  it  would  be  inappropriate  to  use  the  measured  matrix 

conductivity  in  equation  8.  Certainly  the  composite  conductivity  of  the 

1-D  fibre  "F"  composite  (3.0V/  cm  1  K  1  at  300K)  cannot  be  accounted  for 

in  terms  of  the  measured  matrix  conductivity  (0.85W  cm  1  K  1  at  300K) 

and  the  reported  conductivity  of  fibre  "FP  (0.6W  cm  1  K  1  at  300K  e 

The  conductivity  anisotropy  ratios  of  the  1-D  composites  graphitized 

at  2750°C  are  approximately  11-12  at  room  temperature,  decreasing  to 

about  8  at  2300K.  These  are  high  values  and  irrespective  of  whether 

the  bulk  matrix  properties  are  representative  of  the  composite  matrix,  the 

implication  must  be  that  fibre  anisotropy  is  also  high.  It  has  been 
(  22 ) 

reported  that  matrix  conductivity  is  usally  below  both  axial  and 

transverse  fibre  conductivity.  This  is  a  surprising  observation  and  are 

variance  with  our  measured  data. 

Various  models  have  been  put  forward  in  an  attempt  to  explain 

the  transverse  thermal  conductivity  of  a  unidirectional  composite  in 

terms  of  the  thermal  conductivities  of  the  matrix  and  included  fibres. 

(23) 

Some  of  these  have  been  critically  assessed  by  Pilling  et  al  for 
fibre/epoxy  composites.  We  adopt  a  modified  Ohm's  law  approach,  the 
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in-phase  shear  field  analogy,  using  a  model  and  analysis  developed  by 
(24) 

Springer  and  Tsai  .  This  yields  a  generalised  equation  for  the  1-L) 

X 

transverse  thermal  conductivity  X 

J  c 


t 


;d 

_ djr _ 

(2a -6)  +  m/A_  ) 
m  r 


(9) 


o 

where  0  =  f(y),  a  function  relating  fibre  width  at  any  givoiy  and  X  is 

the  fibre  transverse  conductivity.  The  solution  of  Equation  (9)  is  shown 

graphically  in  Figure  66,  for  a  fibre  of  circular  cross-section  and  square 

packing,  as  a  function  of  fibre  volume  fraction  (FVF)  and  X,.  /X  .  Equation 

r  m 

5  assumes  good  contact  between  fibre  and  matrix  but  Figures  34-37 ,43-57  have  shown 
this  assumption  to  be  quite  appropriate  for  the  1-D  composite. 

The  fibre  volume  fraction  of  the  1-D  composites  of  fibre  "A"  is 
50%,  that  for  the  1-D  composite  of  fibre  "F"  is  52%.  Using  values  for 
matrix  conductivity  calculated  from  diffusivity  measurements  and  values 

determined  for  thermal  conductivity  in  a  direction  perpendicular  to  fibre 

X 

orientation  X  c  we  deduce,  using  data  at  350K  and  1000K. 


TABLE  3 


Temperature 

350K 

1000K  1 

spec  type 

X 

X 

CA 

m 

X 

Xf 

r  /X 

m 

X 

X 

c  /X 

m 

X 

Xf 

m 

Fibre  "F"  no  CVD 

0.33 

0.1 

0.38 

0.15 

Fibre  "A"  CVD 

0.12 

<0.01 

0.21 

0.01 

Fibre  "A"  no  CVD 

0.23 

0.03  * 

0.26 

0.05 

X 

X 

c  / 

There  are  indications  that  the  ratio  X^  is  temperature  dependent 
for  all  1-D  composites. 

If  we  now  considered  this  data  in  more  detail  then  three  factors 
can  affect  the  thermal  conductivity  of  the  1-D  composites. 

a)  The  density  of  matrix  graphite  in  the  1-D  composite  is 
higher  than  that  of  the  matrix  graphite  as  supplied. 


b)  The  fibres  within  the  composites  will  have  been  graphitised 
to  a  temperature  of  2750°C. 

c)  Bonding  of  matrix  graphite  to  the  fibre  preform  can  confer 
some  directionality  on  thermal  properties  for  the  composite. 

Matrix  crystallites  in  the  vicinity  of  non  CVD'd  type  "F"  fibres 

.  (17) 

are  oriented  parallel  to  the  fibre  axis. 
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9.  3-D  Composites 

With  the  present  state  of  knowledge  and  lack  of  suitably  relevant 
data  the  modelling  of  C.F.C.C.  can  only  succeed  if  the  problem  of 
quantifying  the  influence  of  processing  can  be  avoided.  This  is  possible 
if  the  3-D  C.F.C.C.  is  considered  to  be  a  two  phase  material  comprising 

a)  Processed  fibre  yarns  representing  the  exial  reinforcement 

b)  Isotropic  matrix  pockets  filling  the  rest  of  the  space  in 
the  unit  cell  geometry. 

9  (  i. )  C .  F  ■  C .  C .  1  s  A  and  C 

The  method  rests  on  the  assumption  that  the  yarns  and  matrix  pockets 

are  arranged  either  in  series  or  in  parallel  and  that  the  resulting  thermal 

model  is  analogous  to  series  and  parallel  connected  electrical  circuits. 

(26) 

Essentially,  it  extends  the  method  used  by  Knappe  and  Martinez-Freire 
to  three  dimensions. 

Along  the  axis  under  consideration,  the  composite  unit  cell  is  sub¬ 
divided  into  four  parallel  conduction  channels  (see  Figure  67;,  it  fs  then 
required  that : 

a)  The  temperature  difference  AT  along  the  heat  flow  direction  is  constant. 

b)  The  total  heat  flow  Q  may  be  divided  into  four  parts  i.e. 

«  =  0,  t  q2  t  Q3  t  q4  ao) 

The  thermal  conductivity  of  each  channel  is  calculated  using  an  Ohm’s 
law  approach.  Thus,  for  the  Z  axis  channels  indicated  in  Figure  67: 


giving 


2aX  X, 
m  lx 

(sX  +  gX  ) 
Ax  m 


where  X is  the  transverse  conductivity  of  the  X  axis  yarn. 
Similarly, 


2aX  X 

Axly _ 

(sX  +  gX  ) 

*  •'  Ay 


(13) 
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2aA  A 
m  lj_ 


3Z  (sA  1  gi  ) 

m  xy 


Lastly, 

X4Z  =  X//Z 


(I'D 


(15) 


where  A.,  is  the  parallel  conductivity  of  the  Z  axis  yarn. 

I !  h 

The  conductivities  of  the  individual  channels  are  then  area-weighted 
and  added  to  obtained  the  total  conductivity  of  the  unit  cell.  For  the  Z 
axis,  the  following  result  is  obtained: 


XlZ^2b-h^r  +  X2Zhr  +  *3Z(2c-r)h  +  X4Z^  2c-r)( 2b-h) 

4bc 


(16) 


And  for  the  X  axis: 


XlXrS  +  X2X^2c“r^S  +  X3X^2c_r^g  +  X4Xrg 
4ac 


(17) 


2bA  A 

where  A^  =  - ...  m-J" V . 


2hxxzxaY 


[hAm  t  (2b-h)^] 


2X 


[  hAlz  +  (2b-h)Axy] 


2bA  A 
Z  m 


3X 


[hA  +  (2b-h) A J 
J.  it  m 


*  X4X  =  X//X 


In  the  real  composite,  these  results  must  represent  only  the  upper 
bound  of  channel  and  unit  cell  conductivities.  Account  must  also  be  taken 
of  the  effect  of  structural  defects  i.e.  split  yarns,  broken  interfaces  and 
porosity,  though  the  latter  may  already  be  adequately  represented  within 
the  1-D  composite  and  bulk  matrix  data. 

For  split  yarns  and  broken  interfaces  perpendicular  to  heat  flow,  the 
lower  bound  to  individual  channel  conductivity  will  be  zero  and  from  this  the 
unit  cell  lower  bound  for  the  particular  defect  type  can  be  found.  A  value 
between  the  bounds  can  then  be  determined  from  the  frequency  With  which  the 
defect  appears  in  the  composite  structure. 

(ii)  C.F.C.C.  'B ' 
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and  perpendicular  to  the  weave  of  a  2-D  CRCC,  using  a  method  due  to 
Bruggeman.  However,  this  approach  uses  the  separately  measured  properties 
of  the  individual  fibres  and  matrix.  Our  investigation  shows  that  these 
properties  are  considerably  affected  by  the  CVD  process  and  fibre  bundle 
densif ication . 

Analyses  of  2-D  woven  composites  have  agreed  the  weave  structure  has 

little  influence  on  the  conductivity  perpendicular  to  the  reinforcement 

plane''  ’  .  Although  the  non  orthogonal  FWPF  construction  increases 

•  ( 23  ) 

conduction  path  length,  this  is  balanced  by  increased  conductivity 

Figure  45  has  shown  that  the  X-Y  weave  is  not  significantly  distorted 
by  the  Z  axis  yarns.  It  seems  appropriate,  therefore,  to  retain  the 
orthogonal  unit  cell  analysis  for  the  Z  axis  conductivity  of  CFCC  'B ' . 

The  case  of  the  parallel  or  in-plane  (X  axis)  conductivity  is  more 

complex.  The  weave  also  results  in  increased  in-plane  conductivity  paths, 

(23)  ... 

but  Kessler  has  estimated  that  the  effect  is  negligible  for  the 

particular  weave  used  in  CFCC  'B'  (eight-harness  satin  weave).  Figures  45 

and  46  have  shown  the  absence  of  clearly  defined  matrix  cross-over  pockets 

in  the  X-Y  weave,  the  yarns  tending  to  fill  the  whole  volume  available. 

The  reduced  incidence  of  the  relatively  high  porosity  matrix  pockets 

may  explain  the  higher  density  of  CFCC  'B'.  Given  that  the  yarn  axial 

conductivity  is  considerably  greater  than  that  of  the  bulk  matrix,  the 

orthogonal  analysis  must  represent  a  lower  bound  to  the  X  axis  conductivity 

of  this  composite.  * 

A  simple  upper  bound  can  be  determined  from  the  assumption  that  the 

heat  flow  paths  within  the  X-Y  weave  are  confined  to  the  X  axis  yarns.  The 

low  conductivity  transverse  paths  through  the  Y  axis  yarns  are  then  effectively 

'short  circuited'  .  Using  the  orthogonal  analysis,  this  may  be  implemented  by 

expanding  the  X  axis  yarn  into  the  volume  normally  occupied  by  the  Y  yarn- 

matrix  conduction  channel.  Although  this  implies  a  reduction  in  yarn  FVF, 

because  of  processing  complications  it  is  not  clear  that  there  must  be  a 

concomitant  reduction  in  yarn  conductivity  such  as  would  be  implied  by 


Equation  'I.  Equation  17  then  becomes: 


_  X2x(2c-r)s  +  X3X(2c-r)c  +  2Xi,xar 
4ac 


33 


10.  Analysis  of  Data 

10 . 1  1-D  C.F.C.C.  Composites 

If  we  follow  the  approach  suggested  in  8.1  to  deduce  thermal 
conductivity  values  for  the  matrix  graphite  within  the  1-D  composite 
then  a  matrix  density  of  'v  1.9  gin  cm  3  would  be  appropriate  for  the  1-D 
fibre  "F"  composite  and  1.67  gm  cm  3  for  the  fibre  "A"  composites.  This 
would  suggest  porosity/tortuosity  factors  of  4.25  and  6.5  respectively 
and  imply  weighting  factors  of  2.1  and  1.4  for  values  of  matrix 
conductivity  to  be  used  in  equations  8  and  9.  Applying  this  to  equation 
8  yields  room  temperatures  values  of  -  4.1  Wcm  3K  1  for  fibre  "F" 

and  Xfn  =3.0W  cm  1  for  fibre  "A"  after  graphitisation  at  2750°C. 

It  is  not  realistic  to  apply  this  approach  to  the  CVD’d  composite  of  fibre 
"  A"  since  we  have  no  value  for  matrix  graphite  graphitised  to  only  2300°C. 

Considering  the  perpendicular  direction  however  and  using  these 
estimated  values  for  matrix  graphite  conductivity  of  1.8  WcnT1  K_1  for 

the  fibre  "F"  composite  and  1.2  Wcm”1  K_1for  the  fibre  "A"  composite 

J. 

yields  values  of  A c  / X m  of  0.15  and  0.167  for  the  two  composites  respectively. 

The  limiting  value  of  A_  /X  <  0.01  would  in  turn  suggest  values  for  X 

r  m  m 

=  1.35  Wcm  1  K  1  in  the  fibre  "F"  composite  and  Am  =  1.0  Wcm"1  K_1  in 
the  Fibre  "A”  composite. 

It  must  firstly  be  asked  if  these  values  are  reasonable.  We  have 
concentrated  on  data  at  low  temperatures  because  the  only  data  available 
for  comparison  is  low  temperature  dat’a  (90-270  K)  obtained  for  carbon 
fibre/epoxy  resin  composites  *  *  .  It  is  reasonable  to  suppose 

that  the  orientation  of  the  "a"  axis  crystallites  of  a  carbon  fibre 

will  be  directly  related  to  its  thermal  history.  The  data  of  Volga  and 

(24)  (28)  . 
Pilling  et  al  clearly  shows  this  to  be  the  case.  Data  by  Volga 

show  a  dramatic  increase  in  thermal  conductivity  of  carbon  fibres 

with  graphitisation  temperature.  The  fibre  thermal  conductivity  A^.11 

at  270  K,  increases  from  1.6  Wcm  1  K  1  after  graphitisation  for  1  hour 


at  2600°C  to  3.0  Wcm  1  K  1  after  1  hour  at  2800°C.  Hence  our  derived 
values  for  X^.11  are  very  much  in  agreement  with  these  observations. 

Likewise  Pilling  et  al  have  determined  the  transverse  thermal 
conductivity  of  UTS  carbon  fibres  manufactured  from  polyacrilonitri te 
precursor  similar  to  fibre  "A"  to  be  0.0G  Wcrn  1  K  1  at  270K.  This  is 
again  in  accord  with  our  conclusion,  that  the  transverse  conductivity 
of  both  fibre  "A"  and  fibre  "F"  is  at  least  one,  and  more  nearly  two 
orders  of  magnitude  lower  than  the  parallel  conductivity. 

The  one  remaining  uncertainty  is  why  the  ratio  X^Am>  using 

assumed  matrix  conductivity  values,  is  less  than  the  limiting  value  of 

0.2  predicted  by  the  model  for  a  F.V.F.  circa  50%.  Hitherto  we  have 

concentrated  on  using  bulk  properties  and  neglected  the  influence  of 

( 15) 

processing.  Stover  et  al  have  shown  that  matrix  crystallites  in 
the  vicinity  of  non  C.V.D.'d  fibre  "F"  composites  are  parallel  oriented 
(P.O.G.)  with  respect  to  the  fibre  axis.  Our  electron  microscopy 
observations  neither  prove  nor  disprove  this  observation.  However  if 
true  this  will  result  in  an  increase  in  effective  fibre  diameter  and 
a  concomitant  reduction  in  the  amount  of  randomly  oriented  matrix  graphite 
in  the  1-D  composite  which  would  adequately  explain  the  low  Xc  /X^  ratio. 

10.2  3-D  C.F.C.C.  Composites 

A  necessary  prerequisite  is  to  establish  the  compositional  similarity 
between  the  1-D  composites  and  the  3-D  yarns.  In  the  case  of  the  1-D 
fibre  a  composite  and  composite  "C"  this  is  patently  inappropriate  since 
the  fibres  in  composite  "C"  are  clearly  of  a  different  type.  Nontheless 
since  we  have  no  accurate  data  for  this  composite  we  will  proceed,  for 
reasons  that  will  be  made  clear  later  on  the  assumption  that  we  can  apply  fibre 
"A'' 1-data  to  composite  "C" 

For  the  general  case  the  yarn  fibre  volume  fraction  is  calculated 
from  the  cross  sectional  area,  the  number  of  filaments/yarn  and  the  affective 
filament  area.  This  yields  a  F.V.F.  of  0.60  -  0.64  ^^which  is  higher  than 
the  F.V.F.  of  either  1-D  composite.  However  the  difference  in  processing  will 
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result  in  a  significant  difference  in  properties  thereby  complicating  any 

attempt  to  model  the  results.  Since  no  experimental  information  is 

available  to  quantify  the  differences  between  CVD  versus  non  CVD  on  fibre 

"F"  composites,  any  attempt  to  determine  the  conductivity  of  a  CVD'd  yarn 

from  the  1-D  composite  data  can  have  doubtful  merit  only.  However,  since 

there  are  no  1-D  composites  with  fibre  "F",  CVD  and  standard  process 

(6) 

available  there  is  no  alternative. Similar  arguments  do  not  now  apply  to 

composite  "C". 

For  fibre  "F"  3-D  composites  it  has  been  shown  that 

a)  Sheath  zones  primarily  of  P.O.G.  can  form  around  the  fibres. 

b)  Where  filaments  are  closer  together  than  about  lg  diameters  the 
space  may  be  entirely  filled  with  C.V.D.  graphite. 

c)  For  distances  greater  than  this  up  to  7  fibre  diameters  the  space 
is  filled  with  transversely  oriented  pitch. 

d)  At  distances  greater  than  approximately  7  diameters  the  orientation 
becomes  random. 

It  is  further  noted  that  where  sheath  zones  predominate  the  matrix  tends 
to  align  the  crystallographic  planes  normal  to  the  filaments.  This  presents 
a  somewhat  idealised  concept  of  the  unit  cell  of  a  3-D  composite  as  shown 
in  Figure  65.  Broken  yarns  are  likely  to  be  matrix  filled. 

Our  observations  on  the  1-D  fibre  composites  indicate  that  the  non 
CVD'd  composites  of  fibre  "A"  and  "F"  show  fairly  dense  composites  with 
little  evidence  of  preferred  orientation  in  the  matrix  graphite.  However 
the  electron  micrographs  of  the  CVD'd  composite  of  fibre  "A"  shows  the  formation 
of  a  distinct  sheath  around  individual  fibres  and  evidence  for  the  formation 
of  transversely  oriented  graphite  in  the  region  between  fibres.  The  density 
of  fibres  is  too  high  to  determine  whether  the  orientation  of  the  matrix  • 
graphite  would  randomise  at  increased  fibre  separation. 

It  therefore  seems  logical  that  the  change  of  orientation  from 
P.O.G.  to  T.O.G.  should  significantly  increase  the  transverse  matrix 
conductivity.  Again  our  arguments  demonstrate  that,  although  the 


thermal  conductivity  of  matrix  graphite  is  higher  than  that  of  the  fibres 
in  the  transverse  direction  the  latter  will  predominate.  Tor  example 
figure  66  shows  that  a  100a  increase  in  matrix  conductivity  will  result  in 
a  35a  increase  in  transverse  conductivity  of  a  composite  with  a  similar 
F.V.F.  A  further  complication  concerns  cracks,  the  density  of  which  is 
higher  in  the  3-D  composites  than  in  the  1-D  composite.  The  majority  of 
these  run  parallel  to  the  fibre  axis.  They  should  have  little  effect 
on  the  transverse  thermal  conductivity  which  is  dominated  by  the  matrix 
but  have  a  significant  affect  on  the  axial  thermal  conductivity. 

On  the  basis  of  the  foregoing  discussion,  the  1-D  composite 
conductivity  data  is  modified  in  the  following  way  for  use  in  the  3-D 
composite  model. 

a)  Axial  conductivity  reduced  by  a  factor  of  0.66 

b)  Transverse  conductivity  increased  by  a  factor  of  2.0. 

No  additional  correction  has  been  made  for  the  difference  in  FVF 

noted  earlier.  The  bulk  matrix  is  used  unmodified  as  its  porosity  is  of  a 

^  0  \ 

similar  order  to  that  of  the  cross-  over  pockets'1 
d)  Modelling  Results 


The  appropriate  unit  cell  constants  are  given  in  Table  4. 


Composite 

Unit  Cell  dimensions  mm 

Type 

a 

b 

c 

g 

h 

r  . 

s 

"A" 

0.42 

0.38 

0.38 

0.42 

0.32 

0.32 

0.42 

"B" 

0.127 

0.625 

0.625 

0.127 

0.71 

0.71 

0.127 

"C"  * 

0.38 

0.43 

0.43 

0.38 

0.34 

0.34 

0.38 

*  Used  as  supplied  by  AFML  although  this  is  known  to  be  incorrect. 

Table  4 

Unite  Cell  Constants 


The  predicted  upper  and  lower  bounds  of  C.F.C.C.  "A"  X  axis  conduct iv 
are  shown  in  Tables  5  and  6  together  with  the  experimental  data  as  defined 


by  the  least  squares  function.  The  rms  errors  are  27 %  and  11%  respectively . 
Clearly  the  lower  bound,  representing  the  case  of  one  broken  yarn/yarn 
interface  per  unit  cell,  provides  the  best  fit  with  the  experimental  data. 

It  was  apparent  from  the  microstructural  investigation  (see  Figures  38  and  39) 
that  the  lower  bound  assumptions  more  closely  reflected  the  real  material. 

The  predicted  upper  and  lower  bounds  of  the  CFCC  "A"  Z  axis  is  shown 
in  Tables  7  and  8.  The  rms  errors  are  12%  and  9%  respectively.  The  results 
are  very  similar  to  those  of  the  X  axis  with  the  greatest  error,  in  the  case 
of  the  lower  bound,  occurring  at  the  extremes  of  the  temperature  range. 

The  predicted  conductivity  of  CFCC  "B"  X  axis  is  shown  in  Tables 
9  and  10.  As  discussed  earlier,  the  upper  and  lower  bounds  of  this  material 
are  different  from  those  of  CFCC  "A".  The  rms  errors  are  10%  and  20%  respect¬ 
ively. 

As  expected  the  lower  bound,  which  represents  the  orthogonal  unit 
cell  geometry,  consistently  underestimates  the  X  axis  conductivity.  The  upper 
bound,  which  though  exaggerating  the  weave  influence,  is  in  much  closer 
agreement  with  the  experimental  data. 

The  calculated  upper  and  lower  bounds  of  CFCC  "B"  Z  axis 
conductivity  are  shown  in  Tables  11  and  12.  The  rms  errors  are  24%  and 
14%  respectively.  The  lower  bound  chosen  is  again  that  of  one  broken  yarn/ 
yarn  interface  per  unit  cell  and  this  seems  reasonable . in  view  of  the  evidence 
provided  by  the  microstructural  investigation  (See  figures  45  and  46). 

For  C.F.C.C.  "C"  we  have  used  in  the  model  data  for  fibre  "A" 

> 

composites  although  it  is  clear  that  the  fibres  in  composite  "C"  care  not 
fibre  "A"  but  more  probably  fibre  "F".  In  Tables  13  and  14  the  calculated 
upper  and  lower  bounds  are  shown  for  the  X  axis.  The  rms  errors  are  6.6% 
and  22.8%  respectively.  The  calculated  conductivities  for  the  Z  axis  are. 
presented  in  Tables  15  and  16  respectively  and  show  r.m.s  errors  of  17.5% 
and  23.8%  respectively. 
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Not  surprisingly  the  modelling  results  for  composite  "C"  are 
poorer  than  those  for  composites  "A"  and  "B".  The  lower  bound  which  is 
arguably  the  most  reliable,  in  view  of  the  microstructural  evidence, 
and  certainly  gives  a  better  prediction  for  composite  "A"  and  the  Z 
axis  of  composite  "B",  consistently  underpredicts  the  thermal  conductivity. 
Clearly  higher  conductivities  particularly  parallel  to  fibre  orientation 
would  be  appropriate.  Comparison  of  data  for  non  CVD'd  1-D  fibre  "F"  comp¬ 
osites  (figures  24-27)  and  fibre  "A"  1-D  composites  (figures  28-31)  show 
the  former  to  be  clearly  higher  by  some  50%  at  350K  decreasing  to  25% 
at  1500K.  If  as  we  believe  composite  "C"  is  also  a  fibre  "F"  composite 
data  based  on  figures  24-27  would  be  more  appropriate  in  the  model  and  lead 
to  better  predictions. 


6. 


Conclusions 


Modelling  of  C.F.C.C.  thermal  properties  can  only  succeed  if 
accurately  quantifiable  information  is  available  on  the  properties  of 
the  individual  phases.  In  this  case  it  has  not  been  possible  to  rigorously 
test  the  present  CFCC  model  because  of  the  unavailability  of  any  suitably 
processed  1-D  composites.  However,  in  spite  of  its  basic  simplicity 
and  the  limitations  of  the  1-D  composite  data,  the  predictions  of  the 
CFCC  model  have  shown  a  reasonable  agreement  with  the  experimental  data. 

For  composites  manufactured  from  fibre  "F"  material,  with  the 
exception  of  the  X  axis  of  composite  "B"  which  is  a  special  case  CFCC 
conductivities  are  reflected  more  closely  by  the  lower  bound  solutions. 

It  is  believed  that  the  microstructural  investigation  has  shown  that  these 
in  turn  reflect  most  accurately  the  structural  state  of  the  real  composites. 
In  the  case  of  composite  "C"  microstructural  evidence  comparing  this  with 
1-D  composites  from  fibre  "A"  leads  to  the  conclusion  that  fibres  are  not 
type  "A"  but  more  likely  fibre  "F".  This  is  reflected  in  modelling  values 
which  consistently  under  predict  and  it  is  suggested  that  better  values 
might  be  obtained  by  using  fibre  "F"  1-D  composite  data. 
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roo.o 

.033 

1.007 

20. R 

650.0 

.MB 

.980 

19. 8 

700. 0 

•  B03 

.951 

JB.7 

750.0 

.788 

.927 

37.6 

roo.o 

.773 

.902 

16.6 

R50.0 

.759 

.877 

15.5 

oon.o 

.745 

•  B54 

14.5 

950.0 

.732 

.831 

13.5 

aooo.o 

.719 

.809 

12.5 

losn.n 

.706 

.788 

11.6 

1100. o 

.594 

.768 

10.7 

1150.0 

.582 

.746 

9.R 

12no.o 

.570 

.730 

9.0 

1?50.0 

.559 

.712 

*.? 

lion.n 

.548 

.696 

7.4 

1150.0 

.537 

.680 

6.7 

3400.0 

.527 

.665 

6.1 

1450.0  - 

.617 

.651 

5.5 

1500.0 

.507 

.638 

5.0 

1550.0 

.598 

.625 

6.5 

1500. 0 

.590 

.613 

4.0 

1550.0 

.591 

.60? 

3.6 

1700.0 

.573 

.59? 

3.2 

1750.0 

.566 

.582 

2.9 

lRoo.n 

.558 

.571 

2.6 

1050.0 

-552 

.564 

2.3 

1900.0 

.545 

.556 

2.1 

1950.0 

.539 

.549 

1.8 

pooo.o 

.533 

.541 

1.5 

2050.0 

•  52B 

.534 

1.2 

2100.0 

.523 

«  5?8 

.9 

2150.0 

.518 

.621 

.6 

2900.0 

.514 

.515 

.2 

2250.0 

.51  1 

.509 

-.3 

2100.0 

.507 

.501 

-.9 

2350.0 

.504 

.496 

-1  .6 

2400.0 

.502 

.490 

-2.4 

2450.0 

.499 

•  4B1 

-3.3 

2500.0 

.498 

.478 

-4.4 

RMS  FRPOR=  12.4 


Table  7 


&  Ul  t*  f  ><J  U  M  7>  >0  ry)  J1  TO 


conluciwjty  hocci.i Wr-«  r>r  <a  5 

l  7  txi c 
LOWER  HOL«ND 


7ri«P£5J6TU°f  (<  J 

ino.n 

isn.o 

4on.n 

450.0 
son.o 
550.0 
6oo.n 
550.0 
7on.o 
75n.n 
nno.o 
850.0 
900.0 
950.0 
inon.o 
lnsn.o 
1 l on.o 
1150.0 
l?nn.o 
J?50.n 
nnn.n 
1150.0 
14.00.0 
14.50.0 
1500.0 
1550.0 
1500.0 
1550.0 
1700. 0 
1750.0 
1500.0 
1850.0 
1900.0 
1950.0 

?non.o 

P050.0 
?l 00.0 
P150.0 
??nn.o 
??50.0 
P100.0 
P150.0 
?4.no. o 
P450.0 
P500.0 


RMS  FR9f>P=  P.9 


r*o.  data 

C*L.  ItTA 

.914 

1.09? 

.916 

1.05? 

.599 

1  .03? 

.58? 

1.004 

.866 

.975 

.549 

.948 

.533 

•  9?1 

.518 

.895 

.503 

.870 

.789 

.845 

.773 

.8?? 

.759 

.799 

.745 

.777 

•  73P 

.756 

.719 

.716 

.706 

.716 

.694 

.697 

.58? 

.680 

.570 

.66? 

•  *59 

.  5£*6 

.548 

.631 

.537 

.616 

•  5?7 

.60? 

.517 

.599 

.507 

.577 

.599 

.565 

.590 

.554 

.581 

.544 

.573 

.534 

.566 

.  5?5 

.559 

.516 

.55? 

.508 

.  545 

.501 

.53? 

.491 

.533 

.487 

.5?9 

.400 

.5?3 

.474 

.518 

.468 

.514 

.46? 

.511 

.456 

.507 

.450 

.504 

.404 

.50? 

.438 

.499 

.43? 

.496 

.4?6 

C  t 


» r  drop 
1*.9 
1^.9 

15.5 

13.5 
IP. 7 
1  1  .5 

10.5 
9.5 

5.4. 

7.3 
6.1 

5.3 

4.3 

3.3 
P.3 

1.4 
.5 

-.1 
-1.1 
-1.9 
-P.5 
-1.1 
-1.9 
-4..  5 
-^.0 
-5.5 

—  5,  0 

—  6  •  4* 
-5.5 
-7.? 
-7.5 
-7.B 
-5.) 
-8.4 
-8.7 
-9.1 
-0.4 
-9.5 

-10.? 
-10.7 
-1  1  .•* 
-11.9 
-1  P.5 
-11.5 
-14.5 


Table  6 


THERMAL  CONDUCT*  Wilt  Hf>OFLlN5  07  C^CC  8 
R  x  axis 
UPPER  BO'»MD 


TEMP£RATURE(KI 

FXP.  DATA 

CAL.  DATA 

*LPPOP 

300.0 

1.3B3 

1.529 

10.5 

350.0 

1.332 

1 .463 

11.4 

<►00.0 

1.263 

1.43B 

12.1 

450.0 

1.237 

1.395 

12. B 

500.0 

1.194 

1.354 

13.4 

550.0 

1.153 

1.313 

13.9 

600.0 

1.115 

1.274 

14.3 

650.0 

1.079 

1.237 

14.6 

700.0 

1.045 

1.201 

14.9 

750.0 

1.014 

1.166 

15.0 

600.0 

.984 

1.132 

15.1 

650.0 

.957 

1.100 

15.0 

900.0 

.931 

1.070 

14.9 

950.0 

.907 

1.040 

14.7 

1000.0 

•  RB4 

1.012 

14.4 

1050.0 

•  B63 

.985 

34.1 

1100.0 

•  644 

.959 

33.7 

1150.0 

.825 

.935 

13.2 

1200.0 

.808 

.911 

12.7 

1250.0 

.793 

.RB9 

12.2 

1300.0 

•  77B 

•  86B 

31.6 

1350.0  _ 

.764  _ 

.848 

10.9 

1400.0 

.752 

.B29 

10.3 

1450.0 

.740 

.811 

9.6 

1500.0 

.729 

.794 

8.9 

1550.0 

.719 

.778 

8.2 

1600.0 

.710 

.763 

7.5 

1650.0 

.701 

.749 

6.8 

1700.0 

.693 

-735 

6.2 

1750.0 

•  6B5 

.723 

5.5 

1600.0 

.678 

.711 

4  .  B 

1650.0 

.671 

.699 

4.2 

1900.0 

.664 

•  6B8 

3.6 

1950.0 

•  65B 

.678 

3.0 

2000.0 

.653 

.  66B 

2.4 

2050.0 

.647 

.659 

1.8 

2100.0 

.642 

.650 

1.2 

2150.0 

.637 

.641 

.7 

2200.0 

.632 

.632 

.1 

2250.0 

.627 

.‘624 

-.4 

2300.0 

.622 

.616 

-1.0 

2350.0 

.617 

.607 

-1.6 

2400.0 

.612 

.599 

-2.2 

2450.0 

.  60B 

.590 

-2.9 

2500.0 

.603 

-5B1 

-3.6 

RMS  ERROR=  10.1 


Table  9 


THtf'.XAl  eOfc/dUCT7*MTY  MODELING  Ur  C *CC  P 

e  *  amc 

LOW£F,  POUND 


TTMP£Rft7  jpe;  IK) 

rxp.  data 

C*l.  DMA 

fcER&OP 

3on.p 

1.763 

1.0H2 

-21.8 

350.0 

1.332 

1.055 

-20.7 

4  00.0 

1.283 

1.029 

-19.8 

450.0 

1.237 

1.003 

-18.9 

500.0 

1.194 

.978 

-18.1 

550.0 

1.153 

.952 

-17.4 

600.0 

1.115 

.927 

-16.8 

650.0 

1.079 

.903 

-16.3 

700.0 

3.045 

.679 

-15.9 

750.0 

1.014 

.856 

-15.6 

flOO.O 

•  9B4 

.833 

-15.4 

850.0 

.957 

.811 

-15.2 

900.0 

.931 

.790 

-15.2 

950.0 

.907 

.769 

-15.2 

1000.0 

•  RB4 

.749 

-15.3 

1050.0 

.863 

.730 

-15.5 

1100.0 

.844 

.711 

-15.7 

115^.0 

.825 

.694 

-16.0 

lPon.o 

.809 

.677 

-16.3 

1250.0 

.793 

.661 

-16.7 

1300.0 

.779 

.645 

-17.1 

1350.0 

.764 

.631 

-17.5 

1400.0 

.752 

.617 

-17.9 

1450.0 

.740 

.604 

•19.4 

1500.0 

.729 

.592 

•18.  p 

1550.0 

.719 

.581 

•19.3 

1600.0 

.710 

.570 

•19.7 

1650.0 

.701 

’  .560 

-20.1 

1700.0 

.693 

.550 

-20.5 

1750.0 

.  6B5 

.542 

-20.9 

1600. 0 

.679 

.533 

-21.3 

1850.0 

.47] 

.526 

-21.7 

1900.0 

.664 

.518 

-22.0 

1950.0 

.659 

.512 

-2  2.3 

2000. 0 

.653 

.505 

-22.6 

2050.0 

.647 

.499 

-22.9 

2) 00. 0 

.64? 

.493 

-23.2 

2150.0 

.637 

,4B7 

-23.4 

■  2?00.0 

.632 

•  4B? 

-23.7 

2250.0 

.627 

.476 

-24.0 

2300.0 

.6?? 

.471 

-?4.3 

2350.0 

.617 

.465 

-24.7 

2400.0 

.612 

.459 

-25.1 

2450.0 

.606 

.452 

-25.6 

2500.0 

.603 

.446 

-26.1 

PMS  ERPnR=  )9.9 


Table  ID 


T Hf KHAL  C0-1'«iCHVJTy  hODFi.  1 N',  (>«■  CFCC 
«  z  amb 

UPPi R  BOUND 


TFmP£RaTUPE  I*) 
300.0 
350.0 
400.0 
450.0 
500.0 
550.0 
600.0 
650.0 
700.0 
750.0 

eoo.o 

850.0 
900.0 
950.0 
1000.0 
1050.0 
1100.0 
1150.0 
1200.0 
1250.0 
1300.0 
1350.0 
1400.0 
1450.0 
1500.0 
1550.0  - 
1600.0 
1650.0 
1700.0 
1750.0 
1BD0.0 
1850.0 
1900.0 
1950.O 
2000.0 
2050.0 
2100.0 
2150.0 
2200.0 
2250.0 
2300.0 
2350.0 
2400.0 

2450.0 

2500.0 


PMS  EPROR=  24.1 


FXP.  data  cal.  data 


.641 

.904 

.633 

.885 

.625 

.866 

•  616 

.847 

•  607 

•  828 

•  597 

.808 

.587 

.789 

.576 

.770 

.565 

.751 

.555 

.732 

.544 

.714 

.533 

.696 

.522 

.678 

.512 

.661 

.502 

.645 

.492 

.629 

.483 

.613 

.474 

•  59B 

.465 

.584 

.457 

.571 

.449 

.558 

.442 

.546 

•  436 

.534 

.430 

.523 

.425 

.513 

.420 

.503 

.416 

.494 

.413 

.486 

.410 

.478 

.407 

.471 

.406 

.465 

.404 

.458 

.404 

.453 

.403 

.447 

.403 

.442 

.404 

.437 

.405 

.433 

.406 

.428 

.407 

.424 

.409 

.419 

.410 

.415 

.412 

.410 

.413 

.405 

.414 

.400 

.415 

.394 

7 ► PPOR 
41.0 

39.7 
38.5 

37.4 

36.4 

35.4  * 

34.5 

33.6 

32.8 
32.0 

31.3 

30.6 

29.8 

29.1 

28.5 

27.8 

27.1 

26.4 

25.6 

24.9 

24.1 

23.3 

22.5 

21.7 

20.8  ' 

19.8 

18.9 
17. R 
16.8 
15.7 

14.6 

13.4 

12.2 

10.9 

9.6 
8.3 
6.9 
5.5 

4.1 

2.7 

1.2 
-.4 

-1.9 

-3.5 

-5.2 


Table  11 


JMLRH&'-  C  OMC  L#c  1  f  v  1  T  r  MODEi.ju'-,  yr  Cf  CC  « 
8  l  A*l< 

ICW£R  POUND 


T£mP£RATUPF  IK) 

rxp.  Data 

CAL.  DATA 

if  PROP 

300.0 

.64) 

.720 

12.4 

350.0 

.633 

.703 

11.0 

400.0 

.625 

.686 

9.7 

450.0 

•  616 

.668 

8.5 

500.0 

.607 

•  651 

7.4 

550.0 

.597  - 

.635 

6.3 

600.0 

.587 

.6)8 

5.4 

650.  P 

.576 

.602 

4.5 

700.0 

.565 

.586 

3.6 

750.0 

.555 

.570 

2.8 

600.0 

.544 

.555 

?.l 

B50.0 

.533 

.540 

1.3 

900.0 

.522 

.526 

.6 

950.0 

.512 

.512 

-.0 

1000.0 

.502 

.499 

~  -.7 

1050. P 

.492 

.485 

-1.4 

1100. P 

.483 

.473 

-2.0 

1150.0 

.474 

.461 

-2.7 

1200.0 

.465 

.449 

-3.3 

3250.0 

.457 

.439 

-4.0 

1300. P 

.449 

.428 

-4.7 

1350.0 

.442 

.  4 1  B 

-5.4 

1400.0 

.436 

.409 

-6.1 

1450.0 

.430 

.400 

-6.9 

1500.-0 

.425 

.392 

-7.7  - 

1550.0 

.420 

.384 

-P.5 

1600.0 

•  41b 

.377 

-9.4 

1650.0 

.413 

.370 

-10.3 

1700.0 

.410 

.364 

-11.2 

1750.0 

.407 

.358 

-12.2 

1  BOO. 0 

•  40b 

.352 

-13.2 

1 850.0 

.404 

.347 

-14.? 

1900.0 

.404 

.342 

-15.3 

1950.0 

.403 

.337 

-16.4 

2000.0 

.403 

.333 

-17.5 

2050.0 

.404 

.329 

-IB. 6 

2)00.0 

.405 

.325 

-19.7 

2150.0 

•  40b 

.321 

-20.9 

2200.0 

.407 

.317 

-22.0 

2250.0 

.409 

•  3)  4 

-23.2 

2300.0 

.410 

.310 

-24.4 

2350.0 

.412 

.306 

-25.6 

2400.0 

.413 

.302 

-26.8 

2450.0 

.414 

.298 

-2B.0 

2500.0 

.415 

.294 

-29.3 

PMS  FPPORs  13.6 


Table  12 


THERMAL 

conductivity 

MODELING  or  crcc 

C 

C  X  AXIS 

UPPER 

BOUND 

TEMPERATURE  <K> 

EXP.  DATA 

CAL.  DATA 

TERROR 

350.0 

.624 

.658 

5.5 

400.0 

.616 

.647 

5.0 

450.0 

.608 

.636 

4.6 

500.0 

.600 

.625 

4.1 

550.0 

.591 

.613 

3.6 

600.0 

.582 

.600 

3.1 

650.0 

.573 

.588 

2.6 

700.0 

.563 

.575 

2.1 

750.0 

.554 

.562 

1.6 

800.0 

.544 

.550 

1.0 

850.0 

.534 

.537 

.5 

900.0 

.525 

.524 

-.1 

950.0 

.515 

.512 

-.6 

1000.0 

.506 

.500 

-1.2 

1050.0 

.497 

.488 

-1.7 

1100.0 

.488 

.477 

-2.3 

1150.0 

.479 

.466 

-2.8 

1200.0 

.471 

.455 

-3.4 

1250.0 

.463 

.445 

-3.9 

1300.0 

.455 

.435 

-4.4 

1350.0 

.448 

.426 

-4.9 

1400.0 

.441 

.417 

-5.3 

1450.0 

.434 

.409 

-5.8 

1500.0 

.428 

.401 

—6.2 

1550.0 

.422 

.394 

-6.6 

1600.0 

.417 

.388 

-7.0 

1650.0 

.412 

.381 

-7.3 

1700.0 

.407 

.376 

—7  •  6 

1750.0 

.403 

.371 

-7.9 

1800.0 

.399 

.366 

-8.2 

1850.0 

.396 

.362 

-8.4 

1900.0 

.393 

.359 

-8.6 

1950.0 

.390 

.356 

-8.8 

2000.0 

.388 

.353 

-9.0 

2050.0 

.386 

.350 

-9.2 

2100.0 

.384 

.348 

-9.3 

2150.0 

.383 

.346 

-9.4 

2200.0 

.381 

.345 

-9.6 

2250.0 

.380 

.343 

-9.7 

2300.0 

.380 

.342 

-9.9 

2350.0 

.379 

.341 

-10.1 

2400.0 

.378 

.339 

-10.3 

2450.0 

.377 

.338 

-10.5 

2500.0 

.377 

.336 

-10.7 

RMS  ERROR*  6.6 


TABLE  13 


thermal  conductivity  modeling  or  crcc  c 

C  X  AXIS 
LOWER  BOUND 


TEMPERATURE (K) 

EXP.  DATA 

CAL.  data 

ftERROR 

350.0 

•  624 

.539 

-13.7 

400.0 

.616 

.530 

-14.1 

450.0 

.608 

.520 

-14.5 

500.0 

.600 

.511 

-14.9 

550.0 

.591 

.501 

-15.3 

600.0 

.582 

.491 

-15.7 

650.0 

.573 

.480 

-16.1 

700.0 

.563 

.470 

-16.5 

750.0 

.554 

.460 

-17.0 

800.0 

.544 

.450 

-17.4 

850.0 

.534 

.439 

-17.8 

900.0 

.525 

.429 

-18.2 

950.0 

.515 

.419 

-18.6 

1000.0 

.506 

.410 

-19.1 

1050.0 

.497 

.400 

-19.5 

1100.0 

.486 

.391 

-19.9 

1150.0 

.479 

.382 

-20.3 

1200.0 

.471 

.373 

-20.7 

1250.0 

.463 

.365 

-21.1 

1100.0 

.455 

.357 

-21.6 

1350.0 

.448 

.349 

-22.0 

1400.0 

.441 

.342 

-22.4 

1450.0 

.434 

.335 

-22.8 

1500.0 

.428 

.329 

-23.1 

1550.0 

.422 

.323 

-23.5 

1600.0 

.417 

.317 

-23.9 

1650.0 

.412 

.312 

-24.3 

1700.0 

.407 

.307 

-24.6 

1750.0 

.403 

.302 

-25.0 

1800.0 

.399 

.298 

-25.3 

1850.0 

.396 

.294 

-25.7 

1900.0 

.393 

.291 

-26.0 

1950.0 

.390 

.287 

-26.3 

2000.0 

.388 

.284 

-26.7 

2050.0 

.386 

.282 

-27.0 

2100.0 

.384 

.279 

-27.3 

2150.0 

.383 

.277 

-27.7 

2200.0 

.381 

.275 

*“f*  «*  •  »i 

2250.0 

.380 

.273 

*■  c  c  .  _• 

2300.0 

.380 

.271 

-28.  ! 

2350.0 

.379 

.269 

-29.1 

2400.0 

.378 

.267 

-29.5 

2450,0 

.377 

.265 

-29.9 

2500.0 

.377 

.262 

-30.4 

RMS  ERROR*  22. B 


TABLE  14 


THERMAL  CONDUCTIVITY  MODELING  OF  CFCC  C 
C  Z  AXIS 
UPPEP  BOUND 


TEMPERATURE <K) 

EXP.  DATA 

CAL.  DATA 

ftERROR 

350.0 

.894 

•  825 

-7.7 

400.0 

.886 

.811 

-8.5 

450.0 

.877 

.797 

-9.2 

500.0 

.868 

.782 

-9.9 

550.0 

.857 

.767 

-10.5 

600.0 

.846 

.752 

-11.2 

650.0 

.834 

.736 

-11.7 

700.0 

.822 

.721 

-12.3 

750.0 

.810 

.706 

-12.8 

800.0 

.797 

.690 

-13.3 

850.0 

.784 

.675 

-13.8 

900.0 

.770 

.660 

-14.3 

950.0 

.757 

•  646 

-14.7 

1000.0 

.743 

.631 

-15.1 

1050.0 

.730 

.617 

-15.5 

1100.0 

.717 

.603 

-15.9 

1150.0 

.704 

.590 

-16.2 

1200.0 

.691 

.577 

-16.5 

1250.0 

.679 

.564 

-16.8 

1300.0 

.666 

.552 

-17.1 

1350.0 

.655 

.541 

-17.4 

1400.0 

.643 

.530 

-17.6 

1450.0 

.632 

.519 

-17.9 

1500.0 

.622 

•  509 

-18.1 

1550.0 

.612 

•  500 

-18.3 

1600.0 

.603 

.491 

-18.5 

1650.0 

.594 

.483 

-18.7 

1700.0 

.586 

.475 

-18.9 

1750.0 

.578 

.468 

-19.1 

1800.0 

.571 

.461 

-19.3 

1850.0 

.565 

.455 

-19.5 

1900.0 

.559 

.449 

-19.7 

1950.0 

.554 

.443 

-19.9 

2000.0 

.549 

.438 

-20.1 

2050.0 

.545 

.434 

-20.4 

2100.0 

.541 

.429 

-20.7 

2150.0 

.538 

.425 

-21.0 

2200.0 

.536 

.421 

-21.4 

2250.0 

.534 

.418 

-21.7 

2300.0 

.532 

.414 

-22.2 

2350.0 

.531 

.411 

-22.7 

2400.0 

.530 

.407 

-23.2 

2450.0 

.530 

•  404 

-23.8 

2500.0 

.529 

•  400 

-24.5 
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TABLE  15 


THERMAL 

CONDUCTIVITY  modeling  of  cfcc 

C  l  AXIS 

LOVER  BOUND 

C 

TEMPERATURE  <K) 

EXP.  DATA 

CAL.  DATA 

*ERROR 

350.0 

.894 

.763 

-14.6 

400.0 

.886 

.750 

-15.3 

450.0 

.877 

.737 

-16.0 

500.0 

.868 

.723 

-16.7 

550.0 

.857 

.709 

-17.3 

600.0  - 

.846 

.695 

-17.9 

650.0 

.834 

.681 

-18.4 

700.0 

.822 

.667 

-18.9 

750.0 

.810 

.653 

-19.4 

800.0 

.797 

.639 

-19.8 

850.0 

.784 

.625 

-20.3 

900.0 

.770 

.611 

-20.7 

950.0 

.757 

.598 

-21.0 

1000.0 

.743 

.584 

-21.4 

1050.0 

.730 

.571 

-21.7 

1100.0 

.717 

.559 

-22.1 

1150.0 

.704 

•  546 

-22.4 

1200.0 

.691 

.535 

-22.7 

1250.0 

.679 

.523 

-22.9 

1300.0 

.666 

.512 

-23.2 

1350.0 

.655 

.501 

-23.4 

1400.0 

.643 

.491 

-23.7 

1450.0 

.632 

.481 

-23.9 

1500.0 

.622 

.472 

-24.1 

1550.0 

.612 

.463 

-24.4 

1600.0 

.603 

.455 

-24.6 

1650.0 

.594 

.447 

-24.8 

1700.0 

.586 

.439 

-25.0 

1750.0 

578 

.432 

-25.2 

1800.0 

.571 

.426 

-25.5 

1850.0 

.565 

.419 

-25.7 

1900.0 

.559 

.414 

-26.0 

1950.0 

.554 

.408 

-26.3 

2000.0 

.549 

.403 

-26.6 

2050.0 

.545 

.398 

-26.9 

2100.0 

.541 

.394 

-27.3 

2150.0 

.538 

.389 

-27.7 

2200.0 

.536 

.385 

-28.1 

2250.0 

.534 

.381 

-28.6 

2300.0 

.532 

.377 

-29.1 

2350.0 

.531 

.373 

-29.7 

2400.0 

.530 

.369 

-30.3 

2450.0 

.530 

.366 

-31.0 

2500.0 

.529 

.362 

-31.7 
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